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Abstract 

We discuss the theory and phenomenology of B ds —> V{—> M\M?)££ decays in 
the presence of neutral-meson mixing. We derive expressions for the time-dependent 
angular distributions for decays into CP eigenstates, and identify the relevant ob¬ 
servables that can be extracted from time-integrated and time-dependent analyses 
with or without tagging, with a focus on the difference between measurements at B- 
factories and hadronic machines. We construct two observables of interest, which we 
call Qg and Qg, and which are theoretically clean at large recoil. We compute these 
two observables in the Standard Model, and show that they have good potential for 
New Physics searches by considering their sensitivity to benchmark New Physics 
scenarios consistent with current b —> sit data. These results apply to decays such 
as B d -f K*{-t K s n°)U, B s -> K S K L )U and B s -> <£(->• K+K~)££. 


1 Introduction 


Rare B decays mediated by flavour-changing neutral currents constitute a unique play¬ 
ground to test the Standard Model (SM) and search for New Physics (NP). Among these, 
processes mediated by the quark-level b —* sii transition have received a great deal of 
attention following a large programme of measurements at R-factories, LHCb and CMS: 
branching ratios, CP asymmetries and angular distributions of B —> 
and B s —> 0/r + /i - [l2| decays, the branching ratio B(B S —> R + R~) (l3-15 , and inclusive 
B —> X s ££ observables 16,17 . Global fits to all b —» sq and b —> s££ data have recently 


uncovered a pattern of tensions between theory and experiment, triggered by the the anal¬ 
ysis of the B —» angular distribution 118-21 , and followed by the measurement 


of the ratio Rk = B(B —> Kfx\x)/B{B —> Kee ), which is consistent with New Physics in 
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R —» K*fifi data and hints at lepton-flavour non-universality 22-28 . In this context, the 


study of new independent b —» observables is of great interest, as a means to gather 

evidence for (or against) these tensions, and to fingerprint the resulting New Physics. 

When dealing with decays of neutral B mesons, experimental observables are affected 
by particle-antiparticle mixing (oscillations), with the decaying meson being either a B 
or a R depending on the time of decay. In the case of flavour-non-specific decays -such as 
decays into CP eigenstates- in which the final state can arise from the decay of both B 
and R mesons, the mixing and decay processes interfere quantum-mechanically, leading to 
interesting phenomenological consequences (for a review see for instance Refs. 29,34 ). In 


particular, new observables arise compared to the case without mixing. These observables 
depend on the experimental set-up (R-factory or hadronic machine), the presence of 
flavour tagging of the decaying R-rneson, and the possibility to perform time-dependent 
measurements (in contrast to the limitation to time-integrated observables). In the case 
of b —y s££ transitions, these effects have been so far taken into account in the untagged 
time-integrated measurements of B s —» 0/r + W 
(see also the discussion in Ref. 


30 
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and B s —>■ /i + /r 
in the case of B s —> VV decays). 


31] at the LHC, 
T ime-dependent 


angular analyses of B^ s —> V££ with tagging are much more challenging experimentally, 
but might be reached at a high-luminosity flavour factory such as Bellc-II 
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In this paper we develop the theoretical framework and study the phenomenological 
advantages of time-dependent B^s —> V££ decays, spelling out the new observables that 
can be accessed, as well as the opportunities for New Physics searches, both at R-factories 
and hadronic machines. While the formalism is valid for any decay of the type B—y 
V(—> M\M- 2 )l£ with M 1 M 2 a CP eigenstate, we identify the following modes of interest: 

• B d -)• K*°{K s n°)£ + £- 

• B s -)■ 0(-> K s K l )£+£- 

• B s -> <j>{-+ K + K~)£+£- 

As a summary of the main points to be discussed below, we shall see that: 

• In the presence of mixing, the time-dependent angular distributions exhibit a new 
type of angular coefficients, hi and Sj, apart from the usual coefficients accessible 


from flavour-specific decays, J* and A, c.f. Eqs. (25),(26). 


Time-integrated CP-averaged rates and CP-asymmetries, as measured at hadronic 


machines, are affected by mixing effects in two ways, c.f. Eqs. (42),(43): 1) The 
terms with Jj± Jj are multiplied by the factors 1/(1 — y 2 ) and 1/(1 +x 2 ) respectively, 
with y = Ar/(2r) and x = A m/T. 2) New contributions proportional to the 
coefficients h t and s t arise. At R-factories only the first type of corrections appear, 
and time-integrated quantities are independent of the coefficients s* and hi, c.f. 


Eqs.(44),(45). 
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• We identify s 8 and s 9 as new observables of interest, which can be extracted most 
conveniently from a time-dependent analysis with flavour tagging. Theoretically 
clean observables can be built from s 8 and s 9 ; two such observables are Q$ and 
Q 9 , which are clean at large hadronic recoil. These observables contain independent 
information, not accessible from flavour-specific decays. Such observables could be 
studied, for instance, at a high-luminosity flavour factory, where a separation be¬ 
tween B and B samples would be possible together with a study of time dependence 
of the decay process. 


• The observables Q 8 and Qg can be predicted in the Standard Model with small 
uncertainties (see Fig. [Tj) . In particular, Qg measures right-handed currents: in the 
case of a b —> s transition, ~ — cos(0 g — 2{3 S ) to a very good precision, with cf) q 
the mixing angle of the B q system. In addition, these observables are very sensitive 
to New Physics scenarios consistent with current b —>• sy and b —> s££ data, such as 
models with Z' bosons with vector and/or axial couplings to fermions. 


The structure of this article is the following. We begin in Section [2] with a discussion 
on time-dependent angular distributions: In Section M we review the basic facts of 
B —> V££ decays without mixing. In Section T2 we address the CP parities associated 
to transversity amplitudes for B —>■ V££ decays into CP eigenstates. In Section 243 we 
derive the expressions for the time-dependent angular distributions, and identify the new 
angular observables h^s) and Si(s ) that arise in the presence of mixing, demonstrating 
in Section 2.4 that S 5 , 6 s, 8 , 9 ( s ) contain independent information not accessible from the 
angular distribution of flavour-specific decays. In Section [3] we discuss in detail the two 
types of observables that can be obtained from the distributions in the presence of mixing: 
time-integrated (Section 3.1) and time-dependent (Section 3.2) observables. We also 
define the observables Q$ and Qg, which are form-factor-independent at large recoil, 
and we provide simplified expressions at the leading order of the effective theory in this 
limit. Standard Model predictions for these observables and New Physics opportunities 
are discussed in Section |4j Finally, we conclude in Section [5] Some details are relegated 
to the appendices. In Appendix |A| we discuss the kinematics of CP-conjugated B —>■ V(— > 
M\MQ££ decays in terms of momentum invariants and the different conventions for the 
kinematic angles that appear in the angular distributions. In Appendix [B] we recall the 
determination of the CP parity for the different transversity amplitudes. In Appendix [C] 
we collect the expressions for the coefficients hi(s ) and Sj(s) in terms of transversity 
amplitudes. 
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2 Time-dependent angular distributions 


2.1 B —> V it decays without mixing 


We first recall a few elements of the analysis of the exclusive b —>• sll decays of the type 
B —> V(—>■ M\M-i)ll. In this subsection, we consider a situation where no mixing occurs, 
and where the M X M 2 state is not (necessarily) a CP-eigenstate. 

This process is described by the usual effective Hamiltonian, with SM operators plus 


(potentially) NP operators with chirality flip, scalar or tensor structure 39 42 


’Hoff — 


4 G f 


V2 L 


\ u CiO't + \ c C 1 O e 1 ~\tY,CiOi 


iei 


( 1 ) 


where X q = V q bV* s and / = {3,4, 5, 6, 8, 7, 7', 9,9', 10,10', S, S', P, P', T, T'}. The opera¬ 
tors (Pi,..,6 and 0 8 are hadronic operators of the type (sTb)(qT'q) and {scji a 'T a P R b)G a 
respectively (431, and contribute to b —y sll processes through a loop coupled to an electro¬ 
magnetic current (via b — > sy* — y sll). These operators are not likely to receive significant 
contributions from NP, as these would show up in non-leptonic B decay amplitude^] The 
operators (P 7 M, 9 «,ioW,sm,p('),xW are given by: 


0 7( ,) = j^ym b [sa^P R{L) b\F^ , 

° 9{,) = . 
° w(,) = , 


0 3 W = J^y[sP R (L)b}[il] , 

o P(l) = ^[s p R{L) b}[i l5 l] , ( 2 ) 

o T « = jFjls^PRmbWa^p^e] , 


with and Pl,r — (1 T 7s)/2- In the SM, and at a scale /x 6 = 0(rrib ), 

the only non-negligible Wilson coefficients regarding the the operators in Eq. (j2| are 
Cf M (/ib) ~ —0.3, Cf M (pb) ~ 4 and Cj’q 1 (//(,) ~ —4 (see Table §; but all might be affected 
by NP. Contributions to B —» Vll from electromagnetic dipole operators 0 7 (i) are (like 
hadronic contributions) of the type b —y sy* —y sll. Contributions from semileptonic 
operators <P 9 (/) 10 (/) si’i p(') t(') are factorizable and their matrix elements can be written as 


(Vll\O sl \B) = (V\T M \B)(ll\T' M \0) , (3) 


where M denotes a collection of Lorentz indices. It is clear that all hadronic, dipole, and 
semileptonic contributions can be recast as decays of the form 


b -»■ y(-» m 1 m 2 )n {-»• ^+r), 


( 4 ) 


1 New Physics contributions at the 
this should have a small impact on b - 


10% level to the operators 0 \, Ch is not excluded 64 . However, 


sll where the effect of semileptonic operators dominates. 
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where N has the quantum numbers of a boson, whose coupling pattern is determined by 
the operators arising in the effective Hamiltonian. In the SM, the structure of 0 7 , O 9, O 10 
shows that N are spin-1 particles, coupling to both left- and right-handed fermions. This 
is in agreement with the presence of 7 * and Z penguin contributions, but it is also able 
to reproduce the contribution from box diagrams involving two W bosons and a neutrino 
((V — A)(V — A) structure in the SM). In an extension of the SM yielding scalar (tensor) 
operators, one should add N bosons with spin 0 (spin 2 respectively). 

We will work under the following assumptions, inspired by the situation in the SM 
and in its most usual extensions 

• CP might be violated in the decay B —>■ VN, but it is conserved in the decay 
N -+£+£-. 


• N can have spin 0 or spin 1, but not spin 2 (no tensor currents in the effective 
Hamiltonian). 

It proves useful to analyse such decays in terms of transversity amplitudes. Let us call 


M mn = 


[US 


(5) 


the hclicity amplitudes for this decay, where m and n denote the polarisations of the 
meson V and the virtual boson N decaying into the dilcpton pair, respectively. 

If N has spin 1, as the initial decaying particle has spin 0, the only combination of 
hclicity amplitudes allowed are ( m,n ) = ( 0 , 0 ), (+ 1 , + 1 ), (— 1 , — 1 ), ( 0 , t), where t denotes 
the timelike polarisation. One can then define the transversity amplitudes (39,,45][49j| 




M. fg+i — M_i_i 

71 


4i 


M + i i+ i + M_ 1,-1 

71 


Aq — M 0 ,o A t — Mqj . ( 6 ) 


The spin-1 N particle couples to the lepton pair either through £ 7 ^ Pi/ or fy M Pi/, and we 
can further separate left- from right-handed components in the amplitudes: Aq, Aq, Ajf, 
H|p A±, A^. On the other hand, due to current conservation and the structure of the 
time-like polarisation e*^{t) oc ( pi+ +Pe-) U , one can see that A t corresponds to a pure axial 
coupling to the lepton pair, vanishing in the massless limit. I 11 the case where N is spin 
0, the only combination of helicity amplitudes allowed is (m,n) = (0,0). The effect of a 
spin -0 particle with a pseudoscalar coupling to leptons can be absorbed into A t , whereas 
a scalar coupling requires a new amplitude, called A s . 


The spin-summed differential decay distribution is given by 45,65 
d 4 r(B -> V(-> M 1 M 2 )e + £~) 


ds dcos 0 M dcos 61 d(f) 


9 

32tt 


J l8 sin 9m + Ju cos 9m + J 2 S sin 9m cos 29 1 


+J- 2 c cos 2 9m cos 29 1 + J 3 sin 2 9 M sin 2 9i cos 20 + J 4 sin 29m sin 29 1 cos 0 


+J 5 sin 29m sin 9i cos 0 + J 6s sin 2 9 M cos 9i + J 6c cos 2 9 M cos 9i 


( 7 ) 


+J 7 sin 29m sin 9i sin 0 + J 8 sin 29m sin 29 1 sin 0 + J 9 sin 2 0 M sin 2 9 L sin 20 , 
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in terms of the invariant mass of the lepton pair s, and three kinematical angles 9g, 9 m, 4> 
(see Appendix [A]). The coefficients of the distribution J, (s) contain interferences of the 
form Re[AjyAy] and Im[AxAy] between the eight transversity amplitudes: 


^oi Aq, An, Ah’', Ay, Ay, A t , As 


( 8 ) 


and are given by 


J Is 


[|Af | 2 + |Af I 2 + |Af | 2 + |Af I 2 ] + ^Re (AyAy* + Af Af*) , 


Jic = \A L 0 \ 2 + |Af | 2 + ^ [|A t | 2 + 2Re(Af Af*)] + ft |A,| 2 , 

J* = f [|^i| 2 + |Af | 2 + |Af | 2 + |Af | 2 ] , J 2c = -ft [|Af | 2 + |Af | 2 ] , 


J* = \& [l^il 2 - I4I 2 + l^l 2 - KH . ^ = ^ft [RfAf Af* + Af Af*)] , 

J 5 = V2 /3, [Re(Af Af - Af Af) - Re(Af A^ + Af *A S ) 

L V s 


J 6s = 2/3, [Re(Af Af - Af Af)] , 


me 


J&c — 4/3, Re(Af A* s + Af A s ), 
V s 


J 7 = f2/3, Im(Af Af* - Af Af*) + -ft Im(Af A£ - Af A s )) 

L \/S 


J 8 = -^/3, 2 [lm(AfAf+ AfAf)] , 


J 9 = ft [lm(Af*Ay + Af Af ] , (9) 


where /3, = fl — 4 m 2 /s. Similar expressions hold for the CP-conjugate decay B —» V(—> 
MiM 2 )P + P-, with angular coefficients J) involving amplitudes denoted by A x , and ob¬ 
tained from the A x by conjugating all weak phases [f} The form of the angular distribution 
for the CP-conjugated decay, however, depends on the way the kinematical variables are 
defined. In the case in which the same conventions are used irrespective of whether the 
decaying meson is a B or a B, we have (see Appendix |a|) : 


dT[B V(-> MiMaf+r] 
ds dcos9g dcos6Rf dcf) 

dT[B -> V{-> M 1 M 2 )£+r] 
ds dcos9g dcos9M d<f> 


'ftft Ji(s)fi(9e, 9 m , <j>) 

i 

^ftC,iJi{s)fi{9e,9M,4>) 

i 


( 10 ) 

( 11 ) 


where fi(9e,9 M ,(f>) are deffiied by Eq. ([7]), and 


Q — 1 for i — Is, lc, 2s, 2c, 3,4, 7 ; Q = — 1 for i = 5, 6s, 6c, 8, 9 . (12) 

for B and B decays, but in agreement with general 


2 This is opposite to the notation used in ref. 
discussions on CP-violation. 


39 
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We stress that this result arises just from the identification of kinematics of CP-conjugate 
decays, and does not rely on any intrinsic CP-parity of the initial or final states involved. 

2.2 CP-parity of final states and decays into CP eigenstates 

The separation into transversity amplitudes not only simplifies the analysis of the inter¬ 
ference pattern, but also provides amplitudes with final states possessing definite CP- 
paritie^j In order to determine the CP-parities associated to the different transversity 
amplitudes we follow the analysis of Ref. [35], where decays of the type B —» MN, with 
M,N unstable particles, are considered. The details of how to apply the results of Ref. [35] 
to the B —$Vlt decays of interest are provided in Appendix [Bj here we briefly summarize 
the main results. 

We consider the decays B — > M\M 2 £ + £~~ and B —>• MiM 2 £ + £~ , such that Mi,M 2 are 
either CP-eigenstates or CP-conjugates, and define the transversity amplitudes: 

A x = A X (B -> MiM 2 £ + £~) , A x = A X (B M x M 2 t£~) , (13) 

where X = LO, RO, L ||, R\\, LX, R_L, t, S. These two sets of amplitudes are related by 

A x = rj x A x (14) 

where rj x are the CP-parities associated to the different transversity amplitudes. We find 
that (see Appendix [B| 

r j X — V for X = LO, L\\, RO, R\\,t ; Vx = ~V f° r X = LX,R±,S, (15) 

where rj = 1 if M\,M 2 are CP conjugates (e.g. K + K~ ), and rj = ~r](Mi)r](M 2 ) if M\,M 2 
are CP eigenstates (e.g. KsKl )• Here rj(M) denotes the intrinsic CP-parity of meson 
M. For the three processes of interest mentioned in the introduction, the combination of 
intrinsic CP-parities leads always to r] = 1. 

At this point we can classify the angular observables J* whether they combine ampli¬ 
tudes with identical or opposite CP-parities, and whether they involve real or imaginary 
parts of interference terms: 

• Real part with identical CP-parities: i = Is, lc, 2s, 2c, 3,4. 

• Real part with opposite CP-parities: i — 5, 6s, 6c. 

• Imaginary part with identical CP-parities: i = 7. 

• Imaginary part with opposite CP-parities: i = 8,9. 

3 We emphasise that the term “CP-parity” makes reference to the final states and not to the amplitudes 
themselves, since the latter involve either a B or a A,which are not CP-eigenstates. 
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We note that the numbers Q defined in Eq. (12) in a different context (identification 
of the kinematics between CP-conjugate decays) corresponds to the product of the CP- 
parities of the amplitudes involved in the interference term Jj. 


We now turn to the case of decays into CP eigenstates: B —> fcp. In this context, it 
is useful to define two different angular coefficients A, Ji which are CP conjugates of Jf 

• the angular coefficients Ji formed by replacing Ax by Ax = Ax(B —>■ fcp ) (with¬ 
out CP-conjugation applied on f C p ), which appear naturally in the study of time 
evolution due to mixing, where both B and B decay into the same final state fcp- 


• the angular coefficients obtained by considering A x = A X (B —» f CP ) (with 
CP-conjugation applied to fcp), which can be obtained from A x by changing the 
sign of all weak phases, and arise naturally when discussing CP violation from the 
theoretical point of view. 


From the discussion above we have Ax = t]xAx, with r)x given in Eq. (15). Plugging 
these amplitudes into the coefficients in Eq. (J9]) , we see that the two types of angular 
coefficients are related through 

Ji = CiJi , (16) 


with Q given in Eq. (12). In addition, in the limit of CP conservation, J, : = Jj. 

Since the final state is not self-tagging, an untagged measurement of the differential 
decay rate (e.g. at LHCb, where the asymmetry production is tiny) yields essentially the 
CP-average 


dV{B f CP ) + dT{B f CP ) 


ds dcosde dcosd M dcj) 


£14 + JiWt, e M , <P) = Ji + CiJiUMt, » m , <t >). 


(17) 

whereas the difference between the two decay rates (which can be measured only through 
flavour-tagging) involves Ji — Ji = Ji — (iJi , 


dT(B f CP ) - dT(B -» f C p) 
ds dcos 9f_ dcos0M d<f> 


</>) = - CiJiWi, Qm, </>). 

i i 

(18) 


We see that the convention chosen in Eqs. @,(0 for flavour-tagging modes allows 
one to treat on the same footing these modes and the modes with final CP-eigenstates, 
since the same combinations of angular coefficients occur in both cases when one considers 
the CP-average or the CP-asymmetry in the decay rate. Let us add that this results from 
a conventional identification between CP-conjugate decays in the case without mixing. 
This freedom is not present in the presence of mixing where both decays result in the 
same final state, which must always be described with the “same” kinematic convention, 
in the sense of a convention that depends only on the final state, without referring to the 
flavour of the decaying B meson (see Appendix [A]) . 







A slightly counter-intuitive consequence is that the CP-asymmetries for J* with i = 
5, 6s, 6c, 8,9 are measured in the CP-averaged rate, and vice-versa. We also note that due 
to the interferences between different decay amplitudes, only some of the Jj — J* differences 
measure CP-violation in specific decay amplitudes (i.e., |A| = |A|, i = Is, lc, 2s, 2c, 3) 
whereas the others measure relative phases between amplitudes (i = 4, 5,6s, 6c, 7, 8,9), 
see Eq. 0 . 


2.3 Angular distributions in the presence of mixing 


In the case of B decays into CP-eigenstates, where the final state can be produced both by 
the decay of B or B mesons, the mixing and decay processes interfere, inducing a further 


time dependence in physical amplitudes (see e.g. Ref. 29,34]). These time-dependent 
amplitudes are given by, 

Ax(t) = Ax(B(t ) —* V(— > fcp ) —> ) = g + (t)Ax + -g-{t)Ax , (19) 

P 


Ax(t) — A x (B(t) —* V{—> fcp)C + P- ) — -g~(t)A x + g + (t)A x , 

q 


( 20 ) 


where the absence of the t argument denotes the amplitudes at t — 0, i.e. in the absence 
of mixing, and we have introduced the usual time-evolution functions 


_ p-^g-r*/ 2 


9+(t) = e 
g-(t) = e- imt e~ Tt/2 


ATt Amt ATt Amt 

cosh-cos- i sinh-sin 


Art 

• sinh —— cos 


2 

Amt 


, Art 

+ i cosh —— sm 


2 

Amt 


( 21 ) 

( 22 ) 


with Am = Mh — M Jj and Ar = — T# (see Ref. |29j). The values of the different 

mixing parameters for the three decays of interest are collected in Table [l} 

In the presence of mixing, the coefficients of the angular distribution also become 


time-dependent, as they depend on the time-dependent amplitudes in Eqs. (19),(20). This 


evolution can be simplified by noting that CP-violation in B q — B q mixing is negligible 
for all practical purpose^, and we will assume \q/p\ = 1, introducing the mixing angle 4 >: 


- = e l4> 
P 


(23) 


This mixing angle is large in the case of the Bj system but tiny for B s , see Table [TJ 
The time-dependent angular coefficients are obtained by replacing time-independent 
amplitudes with time-dependent ones in Eqs. §: 


Ji(t) — Ji[Ax —> /U(l)) . Jiif] — Ji{A_x —> Ax(t)) . 


(24) 


4 The current world-averages are \q/p\B d = 1-0007 ± 0.0009 and \q/p\B d = 1-0038 ± 0.0021 46 . 
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Decay 

V 

0 

sin0 

COS (ft 

Ar 

x = Am/Y 

y = Ar/(2r) 

B d ->■ K*°(K s n°)£+£- 

1 

-2/3 

-0.7 

0.7 

~ 0 

0.77 

0 

B s -)■ ct>{^ K L K s )£ + r 

1 

2/3 s 

0.04 

1 

7^0 

27 

0.06 

B s K+I(-)£+£- 

1 

2/3 s 

0.04 

1 

7^0 

27 

0.06 


Tabic 1: Parameters of the three decays of interest 46 


We consider the combinations J t (t) ± Jj(f) appearing in the sum and difference of time- 
dependent decay rates in Eqs. (Jl7|), (18). From Eqs. (19), ([20]) and (24), we get 


Ji(t ) + Ji(t ) = e Tt (Ji + Ji) cosh(yYt) — hismh(yYt) 


(25) 


Ji{t) 



Ji) cos(aTT) — Si sin(xrt) , 


(26) 


where x = Am/Y , y = AY /(2T), and we have defined a new set of angular coefficients s*, hi 
related to the time-dependent angular distribution. The coefficients Jj, Ji can already be 
determined from flavour-specific decays. The explicit expressions for Sj and hi in terms 
of transversity amplitudes are collected in Appendix [Cj 

Time-dependent angular distributions therefore contain potentially new information 
encoded in the new angular observables Sj and h*. These pieces of information will be 
analysed in the rest of the paper. For the moment a few comments are in order: 


• The coefficients hi are very difficult to extract, since they are associated with 
sinh(yTt) with y very small. In particular, the time dependence of the untagged 
distribution ( Jl7| ) provides essentially no new information. 

• The coefficients s* for i = Is, lc, 2s, 2c, 3,4, 7 are associated with a CP-asymmetry 
in angular coefficients: Ji — Ji. 

• The coefficients S; for i = 5, 6s, 6c, 8,9 are associated with CP-averaged angular 
coefficients: Ji + Ji. 

• The coefficients Sj for i = Is, lc, 2s, 2c, 3,4, 5, 6s, 6c are given by the imaginary part 
of amplitude interferences, s* ~ Im(e ! mjAy), and vanish in the absence of complex 
phases. This is approximately true for B s —> Vii decays in the SM in regions where 
strong phases are small, e.g. in the region s ~ 1 — 6 GeV 2 , and if the NP contribution 
has the same weak phase as the SM. The corresponding coefficients Ji — Ji do not 
vanish, in general. 

• The coefficient s 7 vanishes in the absence of phases in the amplitudes, while the 
combination J 7 — J 7 vanishes in the absence of CP violation in decay. Both are 
therefore very small in the SM, and also if the NP amplitudes have approximately 
the same phase as the SM. 


10 












In the same conditions as above (no complex phases), the coefficients (J t + Ji)i= 8i9 
vanish, while s 8j9 do not. 


It seems therefore that the most promising observables in this context are s 8j9 , which 
could be large and can be extracted from the time evolution of 


Js(t) — J 8 (t) ~ — s 8 e r< sin(a;rf) , J 9 (t) — J 9 (t ) ~ — s 9 e lt sm(xTt) . (27) 


The coefficients s 8 and s 9 have the following expressions (see Appendix [C|: 
1 


s 8 = ~n=Pt Re 


e^A^A? + A*Af) + e-^{A^A L * + A$Af) 


s 9 = $ Re e**(Af:A^* + AKAf) + e-^(Af\A^ + Aff Af) 


(28) 

(29) 


We have checked by direct calculation that indeed the coefficients s* with i ^ 8,9 are tiny 
in the SM, and that they do not get significant enhancement from NP contributions if 
new sources of CP violation are not large. 

We emphasise that the measurement of the coefficients s 8)9 is challenging from the 
experimental point of view, since the study of ■J l {t) — A(f) requires 1) flavour tagging of 
the original sample to separate B and B at t — 0, 2) the use of appropriate foldings to 
extract the corresponding angular contributions, identical to the ones used to extract J 8 
and J 9 12,31, and 3) a time-dependent analysis to isolate the sin(xrt) coefficients. 


2.4 Symmetries of the distribution 


Having identified a few new observables accessible from the time-dependent angular distri¬ 
butions, it remains to be seen if they are truly independent from the observables that can 
be extracted from angular distributions of flavour-specific decays. The information that 
can be obtained from the angular distributions depends on the number of independent 
combinations of interference terms Ax Ay in the angular coefficients. A systematic for¬ 
malism to determine which combinations can be accessed from the angular distributions 
alone is the “symmetry formalism” developed in Refs. 40,47 .^] 

In the approximation of massless leptons, and neglecting scalar and tensor operators, 
the angular distributions of flavour-specific decays contain a unitary symmetry, given by 


the transformation 40 


rii = 


A L 


—> Uru 


(30) 


with U an arbitrary unitary 2x2 matrix, and {co, <7||, cr_i_} = {1,1, —1}. Under this 
group of transformations, J, —> J t . This means that from flavour-specific decays, only 
those combinations of terms Ax Ay that remain invariant under this transformation can 

5 See also Ref. 48 for an application to S- and P-wave components in B —► (Ahr)/z/i. 
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be accessed. This approach is useful to eliminate redundancies among observables built 


from the angular coefficients J* 40,47 


We now identify the transformation properties of the coefficients s* -neglecting weak 
phases for simplicity. We note that, under the unitary transformation: 


R e[A?Af* ± AfAf] —► [1 - (1 t ^)A 2 ] Re[AfAf ± AfAf*] 

+(1 =f an) Re [a^iA^Aj + aj-qAfAf] , 


(31) 


Im[AfAf* ± AfAf] —► [1 - (1 ± aij) A 2 ] Im [AfAf ± AfAf] 

“(I ± °a) Imfr/Af Af - a^AfAf] , (32) 

where oy,- = a t aj, A 2 = 1 — |f/n| 2 , rj = U\\U* 2 and i,j = 0, ||, _L. Non-trivial transforma¬ 
tions involve only Re [AfAf ± AfAf] with a tJ = =pl, or else Im[AfAf ± AfAf] with 
atj = ±1. From the explicit expressions given in Appendix [Cj we see that (neglecting 


lepton mass terms and weak phases in the amplitudes): 

Sis, 2 s ~ sin0 ■ Re [AfAf + AfAf] — (|| —Y _L) , (33) 

sic, 2 C ~ sin 0 • Re [Af.Af : * + AfAf] , (34) 

s 3 ~ sin0 • R e[AjAf* + AfAf ] + (|| —>• _L) , (35) 

s 4 ~ sin0 • Re [AfAf* + AfAf*] , (36) 

s 5 ~ cos 0 • Im [Af Af* — Af Af] , (37) 

s 6s ~ cos0 • Im[Af Af — Af Af] , (38) 

s 7 ~ sin 0 • Im [Af Af* - Af Af*] , (39) 

s 8 ~ cos 0 • Re [AfAf* + Af Af] , (40) 

Sg ~ cos 0 • Re [Af Af* + Af A? ] . (41) 


Therefore the only coefficients s t that (in this approximation) do not remain invariant 
are S5, Sq s , s 8 and Sg, which contain additional information not accessible from the usual 
angular distributions of flavour-specific decays such as Bd —y K*(—> K + Among 
these coefficients, we have seen that s 8j g are particularly promising; now we see that they 
are independent from the coefficients J t . 
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3 Observables 


The expressions in Eqs. (25), (26) for the coefficients of the time-dependent distributions, 


show that additional structures arise in the presence of neutral-meson mixing. In this 
context, two different quantities might be considered: time-integrated observables, or 
observables related to the time dependence. In this section we discuss the two possibilities. 


3.1 Time-integrated observables 

As discussed in Refs. 32 ,[34], time integration should be performed differently in the con¬ 


text of hadronic machines and R-factories. The time-dependent expressions in Eqs. (25) 


and (26) are written in the case of tagging at a hadronic machine, assuming that the two 


5-quarks have been produced incoherently, with t G [0, oo). In the case of a coherent BB 
pair produced at a R-factory, one must replace exp(— Ft) by exp(—T|t|) and integrate 
over t G (—oo, oo) J34j|. Interestingly, the integrated versions of CP-violating interference 
terms are different in both settings, and the measurement at hadronic machines involves 
an additional term compared to the R-factory case: 


(Ji + Ji) Hadronic p 


y 


{Ji - Ji) Hadronic 

{Ji + Ji) B—factory 
{Ji - Ji) B—factory 


1 + X 2 
2 1 


X (Jj + Ji) 


X (Ji Ji) 


y 


X 


1 + X 2 


x hi 


X Si 


n -y 2 
2 1 
fl + x 2 


[Ji + Ji] , 
[Ji - l] ■ 


(42) 

(43) 

(44) 

(45) 


Making contact with experimental measurements requires to consider the total time- 
integrated decay rate. The time-dependent rate is given by 


dT f , 

df~J dt 


~{^Jls{t) + J\c{t )) — -(2J 2s {t) + J 2 c{t)) 


which after time-integration becomes 

dT \ _ 1 

dq^/ ~ T (1 — y 2 ) 


C) , 


{!) 


Hadronic 


2(«/i s “I - Ju y s ) ~\~ (Jie ~\~ J\c y h lc ) 

2 («/2 s “1“ J2s y ^2s) H - («^2c “1“ J%c y ^2c) 


(46) 


(47) 


(48) 
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(X)b— factory — (X)Hadronic (h — 0) , 


( 49 ) 


where X is the usual normalisation considered in analyses of the angular coefficients. The 
normalised time-integrated angular coefficients at hadronic machines or 5-factories are 
therefore: 


(E*) 


Hadronic 


<Si) 


B—factory 


(A,) 


Hadronic 


(A.) 


B—factory 


{Ji ~ I - Ji )Hadronic _ {Ji Ji) V ^ ^li 

(dY /dq ^ 1 } Hadronic \E) Hadronic 


{Ji + Ji) B—factory 
{dT / dq^)j $—factory 


= (E 4 ) 


Hadronic 


(h = 0) 


(Ji ~ Ji) 

(< dT/dq 2 ) 


Hadronic 


Hadronic 


1 - y 2 (Ji, - Ji) - x x Si 

1 -f- X ^ (X) Hadronic 


- Ji) B—factory 
( JT / d<r)B-factory 


(Aj)Hadronic.(^ ^ 0) . 


(50) 

(51) 

(52) 

(53) 


We see that the interpretation of the time-integrated measurements (£j) from dT(B —> 
fcp££) + XT(5 —> fcp££) is straightforward in terms of the angular coefficients at t — 0. 
Even in the B s case, the smallness of y means that h t will have only a very limited impact 
on the discussion. The time-integrated terms (A*) from dT(B —> f C p££) — dY(B —> fcp ££) 
are subject to two different effects, in particular for B s where x is large: 


(a) they receive contributions proportional to x and y with a different combination of 
interference terms (in the case of a measurement at a hadronic machine), 


(b) they are suppressed (in all experimental set-ups) by a factor (1 — y 2 )/( 1 + x 2 ). 

The discussion above applies in particular to the measurement of B s —y <$>(—» K + K~)££ 


as performed at LHCb 12 . Since this is not a self-tagging mode, and assuming that there 


is an equal production of B s and B s , what is measured is dT(B s —> 0(—» K + K~)££) + 
dT(B s —> <f>(—> K + K~)££), so these measurements have access to the following combina¬ 
tions: 


(Ji T Jj)Hadronic Rr i Is, lc, 2s, 2c, 3, 4, 7, 

(Ji Jj)Hadronic fo^ ^ 5, 6s, 6c, 8, 9. (54) 


The time-integrated observables (S 6 ) H adronic an d (Eg) Hadronic have already been measured 
(under the name of A e and Ag) in Ref. 12 , and are indeed measuring CP-violation. In the 
context of the extraction of $$ and Sg at hadronic machines, one expects them to dominate 
(Ag 9 )Hadronic; especially in the case of B s decays where x enhances their contribution with 
respect to the (J* — J,) term. However, they are overall suppressed by a factor ~ 1/x, 
which in the case of B s decays is quite effective (1/x ~ 0.04). In addition, the necessity 
to perform initial flavour tagging makes these measurements very difficult at hadronic 
machines. 
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We see therefore that (Ej) contain essentially the same information as (Ji+Ji), whereas 
(A i) have a potentially richer interpretation, but are suppressed and thus probably diffi¬ 
cult to extract experimentally. In the following section we will see that time-dependent 
observables do lead to more interesting opportunities. 


3.2 Time-dependent “optimised” observables with tagging 

From the discussion in Section 


is clear that a full tagged time-dependent angular 
analysis of the decay B —)■ V[—> (MiM 2 )cp]^ provides a measurement of the angular 
observables J,(s), Jj(s), Sj(s) and h;(s), i.e. Eqs. (17), (18), (25), (26). We have also seen 


that the coefficients hi are the sinh(yrt) coefficient of J*(f) + Jj(t), whose effect remains 
negligible for t < 10 Ts q , constituting a rather difficult measurement. 

From the theoretical point of view, these observables are quadratic in hadronic form 


factors (see Section 4.1). For instance, 


Si{s) ~ A x Ay 


7 b^ V( 


■F 


B->V i 
Y \ 


s) > 


(55) 


and similarly for Jj(s), J ? ;(s), hi(s). Here F xy represent (schematically) hadronic B —* V 
form factors related to the amplitudes A X y. These form factors constitute a major 
source of uncertainty in the theoretical predictions for the observables. This problem 
is usually tamed by defining a class of special observables with reduced sensitivity to 

These “optimised” observables can be constructed 

and 

at low recoil (s ~ m 2 B ) (50,54 , where the use of effective field theories (SCET (59, 


form-factor uncertainties 40,49-54 


and 

HQET 66,671 respectively) ensures a complete cancellation of form factors at the leading 
order in the respective expansion^ 

In the following, we focus on the large-recoil region for definiteness. We consider the 
following optimised versions of the observables s 8i9 : 


Q 8 — 


S 8 


— 2 (J 2 c + J 2 c)[‘ 2 (J 2 s + J' 2 s ) — (J3 + <73)] 


s 9 


(56) 


Q9 — 

2 ( J2s + J2s) 

There are other possible normalizations for s 8 that are also optimised at large recoil: 

Q 8 + = 


(57) 


— 2 (J 2c + J 2 c)[ 2 (J 2 s + J2s) + (J3 + <73)] 


(58) 


6 In the following and throughout the paper we use the term “large-recoil limit’’ to denote the following 
approximation valid in the region s m 2 B : leading order in a s and leading power in the SCET expansion. 
This is, of course, a slight abuse of language. 
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(59) 


Qs = 


S8 


— 2(J2c + ^2c)[2(^2s + ^2s)] 


The observable Q$ has the particularity of being also optimised at low recoil. However, 
we find that both Qf and Qg are slightly less sensitive to NP than Qg . While it might be 
worthwhile to study these observables further, we will focus here on Qg for illustration, 
noting that its properties do not differ much from those of Qf, Qg. 

Concerning Qg, other possible normalizations involve Jq s or Jg, both of which lead to 
observables that are optimised also at low recoil. We do not consider these possibilities 
any further as the denominators contain zeroes within the kinematical region of interest. 


It is useful to consider these observables in the large-recoil limit (see e.g. Refs. 40,41]), 
where the expressions simplify considerably, the cancellation of form factors is exact, and 
the dependence on the Wilson coefficients is apparent. We find: 


Qs — 


C7 (^7 + Cg ) 


2 |C 7 Iy/(C i0 ) 2 + (2 C 7 + Cg ) 2 


(^ 10 ^ 10^7 + (2C 7 +C 9 X CfCg +C 7 Cg)) s 
4Cf|c 7 -|y(cr 0 ) 2 + (2c 7 - + c s -) 2 K 


+ ... 


COS <f) q 


(60) 


Q9 ~ 


2 CfCf 


[(C7 ) 2 _ (C7 ) 2 ] (C7 Cg ~ CyCg) S 


l (c 7 -r + (c+y 


1(C 7 -) 2 + (C 7 + ) 2 ] 2 


m 


+ 


B 


cos <p q (61) 


where we have assumed real C,;, and used the notation Cf = Ci ± C l >. In the case of the 
b — > s processes at hand, we have (j) q = (j) q — 2/3 s , with (f) q the mixing angle in the B q 
system. We note that if C = 0 (that is, Cf = C~), on has Qg = — cos <j) q , so that the 
value of (Q g + cos (j) q ) is a measurement of right-handed currents. 

I 11 the following section we give Standard Model predictions for these observables and 
study briefly their sensitivity to New Physics. 


4 Numerical Analysis 

4.1 Standard Model 


The systematic formalism to B —y Vii decays at large hadronic recoil to NLO in QCD- 


factorisation has been presented in Ref. 55 and is by now quite standard. In our analysis 


we follow closely the procedure of Refs. 54,58 to which we refer the reader for further 


details. The different transversity amplitudes can be written as: 


^_l^( s ) — Ai i (C9 + l^( s )(l + Vpc) + A utY u (s )) 


~F C+ 


V(s) 

M + 


m 


+ 2 m 1 


. (62) 
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CiQ iib) C 2 (nb) C^{nb) C^Hb) Cb(nb) C^(fib) C^ s (/j,b) Cg S (nb) Cg(fib) Cio(nb) 

-0.2632 1.0111 -0.0055 -0.0806 0.0004 0.0009 -0.2923 -0.1663 4.0749 -4.3085 


Tabic 2: Wilson coefficients in the Standard Model at NNLO at the scale fib = 4.8 GeV. 


' R (s) — -A/j| (C g + lf(s)(l + t^pc) + A u tW(^)) TC 10 


( 2 )> 


^o ’ K ( s ) — K (C g + y t (s)(l + r^pc) + A utY u (s)) =F C 10 

2 nib 


Ai(s) 2 m b 

M + m + s 2 r ’ 

^ 12 ( 3 ) 

M + m 


(M 2 + 3 m 2 — s)T 2 


A 


% 


MS) ~ % K" + 2 me 0 ’’ 


2 M 2 — m 2 3 


A)(s) j Aj(s) — — J\f± C s A 0 (s) , 


where M = ms q and m = my, and: 

• The normalizations are given by 

A/1 = V2\N , jVJ| = s/2(M 2 - m 2 )N , Af 0 = -N/(2m^s) , 
with A = M 4 + m 4 + s 2 — 2 (M 2 m 2 + M 2 s + m 2 s), fit = \/l — 4m 2 /s, and 


N(B) = v; b v ts 


G F a 2 s A 1 / 2 ^ 
3 • 2 10 7t 5 M 3 


1/2 


, 7V(R) = E«y t 


£,S 


G 2 F a 2 s\ l / 2 f3 e 

3 • 2 10 7t 5 M 3 


1/2 


(63) 


(64) 

(65) 


( 66 ) 


(67) 


Y t (s) and l„(s) are the 1-loop contributions from 4-quark operators to the photon 
penguin with the structure sometimes combined with Cg into C 9e fj(s). Y t (s) 

denotes the contribution proportional to V t bV t *, and can be found in Eq.(10) of 
Ref. 55 . Y u (s) denotes the CKM-suppressed contribution, which is multiplied by 
the prefactor A ut (B) = V ub V* s /V t bV bs or A ut (B) = V* b V us /V* b V ts . This function can 
be found in Eq. (A.3) of Ref. [56]. 


• The functions 71 encode contributions from dipole operators C 7 ( 0 , and the rest of 
the hadronic contributions not contained in Y t ,Y u : 


V(«) 

= Cf«T 1 { s )+T±(8){l+$!) 

(68) 

r 2 -(s) 

— C 7eS T 2 (s) + m 2 71(s)(l + rj PC ) 

(69) 

%-(s) 

— C 7e ffT 3 (s) + 71 (s) + 7jj (s)(l + rjic) 

(70) 


The quantities 71 11 represent factorizable and non-factorizable hadronic contribu¬ 
tions in QCD-factorisation and can be extracted from the formulae in Section 2 of 


17 
































Ref. [55], They depend on distribution amplitudes and on two “soft” form factors 
£_l(s), £||(s). The “effective” coefficients C± eS include contributions from 4-quark op¬ 
erators with b and s-quark loops with the structure s[<$, ^/Q\b (see e.g. Refs. 44,571). 


The parameters rjpl in Eqs. (]62|)-(|64|) and (|68|)-(|70|) parametrize non-factorizable 


0(A/rrib) power corrections absent in the current QCD-factorisation calculation. 
Following Ref. 


58 , we write 


Vpc = r?e 


= r“e^“ + r^is/M 2 ) + r^(s/M 


2\2 


(71) 


and take r“’ 6,c = 0 as our central value, varying the parameters r*, <f>i within the 


ranges r“’ ,c G [0,0.1] and </>“’ ' c G [—7r,7r] in the error analysis. This corresponds 
to a contribution from non-factorizable 0(A/rrib) corrections of (9(10%) with an 
arbitrary phase. 

The functions V(s), %o(s), Ai(s), A 2 (s), Xi(s), T 2 (s), T 3 (s) represent the seven inde¬ 
pendent B —» V QCD form factors (see e.g. Refs. (59, 60 ), with the combination 


Ai 2 (s) = (M 2 — m 2 — s)(M + m) 2 Ai(s) — AA 2 (s) 


(72) 


entering Aq' r . Following Ref. 54 , we use the large-recoil symmetry relations 


59 


60 


to express these form factors in terms of £_l(s) and £||(s), defined in the “scheme 1” 


of Ref. 58 , including factorizable power corrections. At a second stage, these are 


themselves expressed in terms of R(s), Ai(s), A 2 (s), which are taken from the light- 
cone sum rule calculation of Ref. | 6 l], both for B —» K* and B s —> 0 transitions. 


The Wilson coefficients Cf are defined as: Cf = C, ± C t i. The Standard Model values 
for these coefficients are collected in Table [2j computed at a renormalisation scale 
fib = 4.8 GeV. As has become customary in analyses of B —» Vii decays 18-25 , we 


use NNLO Wilson coefficients, keeping in mind that the NNLO scheme and scale 
ambiguity can only be eliminated by including the (currently unknown) NNLO 
matrix elements. In this context, in the error analysis we consider a variation of the 
renormalisation scale /i G [/ift/2, 2fib]. In addition, we have Cf, M = (m s /mb)Cy M . 


Following this procedure, we compute central values and errors (by means of a flat 
scan over all parameters) for the observables Q$ and Q 9 in the SM. We compute these 
observables differentially in s, keeping in mind that a proper comparison with data would 
require an integration over bin ranges (see for example the discussion in Ref. (531). 

Our SM results for the observables Q$ and Qg in the low-s region are shown in Fig. [l} 
We show both cases: B s (t) —> (f>(—> KK)/i + fi - and B d (t) — > /!*(—)• A% 7 r°)/i + /i _ , noting 
that the results are very similar. We see that indeed the observable Qg ~ — cos (f) q in the 
whole region (with 4>d — —2 (3 and (j) s = 0 in the SM), while Q § features a distinctive 
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Figure 1: SM prediction for the observables Q$ and Qg in the case of B s {t) —>■ 0 (—* 
K + K~)n + n~ (upper row) and B d (t ) —>• K*(—t A’s7r)/i + /i _ (lower row), in the large-recoil 
region, including error estimates from all sources. See the text for details. 


shape with a zero at so — 2 GeV 2 . This is located at the same position as the zero of sg, 
which can be expressed solely in terms of Wilson coefficients taking the large-recoil limit, 


so 


-2Cy~(2Cy T Cg ) sm — 2 C’j{ 2 C'j T Cg) 


m 


B 


^10^10 ffi (2^7 + Cg )Cg Cf Q + (2C7 + Cg)Cg 


(73) 


The position of this zero measures a different ratio of Wilson coefficients compared to the 
zero of other observables, such as /Ifb or P 2 40 . We stress that the bands in Fig. [I] include 
all sources of error including parametric and form-factor uncertainties, as well as our 
estimates of power corrections, exhibiting the theoretical accuracy for these observables 
in the Standard Model. 


4.2 New Physics 

We now study the sensitivity of the observables Q$ and Qg to different models of New 
Physics. We start with a general scan of (real) New Physics contributions to Wilson 
coefficients compatible with all current constraints from rare P-decays, in order to assess 
the NP reach of the new observables. For that purpose, we write 

C t = Cf M +Cr, (74) 
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Figure 2: NP reach of the observables Q 8 and Qg in the large-recoil region. See the text 
for details. 


and consider the 3cr ranges for the NP contributions Cf p (at the scale Hb = O(mb)) that 
were obtained in the global fit to b —> sy and b —> sit data of Ref. fl8]: 


C, 


NP 

7 

NP 


-0.08,0.03) , 


cr o p e(-2.0,3.0), 


C 9 P e (-2.1, -0.2) , 

Cy, p e (-0.14,0.10) , C 9 T e (-1.2,1.8) , e (-1.4,1.2) . 

The result of this scan is shown in Fig. [2} We consider separately three scenarios: 


(75) 


• LHC (Left-Handed Currents) scenario: NP contributions to C 7 ,Cg,Cio only. This 
corresponds to the orange regions in Fig. [2j delimited by dashed lines (along the 
line Qg = —1 on the right-hand plot). 

• RHC (Right-Handed Currents) scenario: NP contributions to Cr,Cgi,Cio/ only. This 
corresponds to the red regions in Fig. [2j delimited by dotted lines. 

• General NP scenario: NP contributions to all six coefficients C 7 w, C 9 o ), C 10 m . This 
corresponds to the regions in green in Fig. [2| with solid borders. 

We also show the SM predictions for comparison (blue bands in Fig. |2j with ~ — 1 
and ~ —0.7 for the B s and Bd cases respectively). We see that NP can indeed have 


20 






























Figure 3: NP benchmarks for the observables Q s and Qg in the large-recoil region. See 
the text for details. Benchmarks A (blue), B (red) and D (dashed) are hardly visible. 


a large impact on , Qg. As discussed in Section 
from — —cos(J) q requires right-handed currents. 

We finish our exploratory NP analysis by studying a few motivated benchmark NP 
scenarios: 


3.2, any significant deviation of Qg 


A. 


Best fit point in the C 7 — Cg scenario of Ref. 



q p = -0.02, c 9 np = -1.6. 


B. Best fit point in the Cg — Cg' scenario of Ref. 



C 9 np = -1.28, C 9 n , p = 0.47 . 


C. Z'-motivated C 9 (/>, C 10 co scenarios (see e.g. Refs. 



ci. c 9 np = -cr 0 p = -1 

P9 pNP _ _pNP _ 1 

L^g/ - ^10 / — ^ 

pq pNP _ pNP _ _pNP _ _pNP 

— ^9' — u 10 — u 10' 


-1 
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P/I r>NP pNP r>NP pNP i 

Ug — Ug, — U 1Q U 10' _ 1 


D. Best fit point in the general fit of Ref. 18 


c* p = -o.o2, c 9 np = -i.3, cf p = o.3, cyr = -0.01 , C 9 n , p = 0.3 , CJ£ = 0 . 

Scenarios C.l and C.2 arise also respectively in singlet/triplet and doublet leptoquark 
models motivated by recent data on the ratio B(B —> —> Kee ) (see Ref. [22]). 

The predictions for the observables Q$ and Qg within the benchmark scenarios are 
shown in Fig. [3j together with the SM prediction. We see that, among the considered 
scenarios, the only ones leading to significant deviations with respect to the SM are 
scenarios C (corresponding to large NP contributions to Cgt and Cio')> while scenarios A, B 
and D are very close to the blue band (corresponding to the SM prediction). As discussed 
before, scenario Cl has no impact on Q 9 as it has no right-handed currents. Therefore, 
measurements of these observables compatible with the SM would give support to the 


best fit points obtained in the global fits of Ref. 18 25 that we have considered, with the 

As 
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potential to exclude the scenarios with Cg ( p , C^ p (/) such as the one discussed in Ref. 
an alternative viewpoint, these observables could test the latter scenarios, and provide an 
alternative confirmation if more accurate measurements for time-integrated observables 
happened to confirm any of them. 


5 Conclusions 


Decays of the type B d , s V(—> MiM 2 )£ + £~ mediated by the underlying flavour-changing 
neutral current process b —y sli are of great phenomenological interest for two reasons: 
they lead to a vast set of independent experimental observables, and they exhibit a remark¬ 
able sensitivity to New Physics. The decay mode B d —y K*°(—y K~7t + )/j, + /i~ has been 
the first one to be carefully scrutinized, both experimentally 


11,68] and theoreti¬ 
cally 18 21 {69 -111 , and first angular analyses of the decays B s —* (j>{—> K + K~)fi + fi~ jl21 
and B d —> K*°(—> K~n + )e + e~ 112 have been already performed. 

In the case where M X M 2 is a CP eigenstate (such as B d —> , B s — * 

<t>(— > K + K~)h + /j,~ or B s —> 0(— » KsKl)/j , + neutral R-meson mixing interferes with 
the decay, leading to interesting differences with respect to flavour-specific processes where 
mixing plays no role (such as B d —» K*°(—> K~ 7r + ) ji + fi~) . In this paper we have studied 
the effects induced by neutral-meson mixing for the analysis of exclusive B —^ V£b decays, 
spelling out the theoretical formalism and analysing its phenomenological consequences. 

As a first observation, the angular distributions become time-dependent, with addi¬ 
tional structures compared to the case without mixing. These structures are the new 


angular coefficients hi and s*, defined in Eqs. (25),(26) and given explicitly in terms of 


the different amplitudes in Appendix 0 Two types of observables can then be defined 
for these modes: time-integrated and time-dependent observables. The first type depend 
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on the experimental set-up (5-factory or hadronic machine) and differ from the corre¬ 
sponding observables in decays without mixing by multiplicative factors depending on the 
mixing parameters x and y. In addition, the expressions for time-integrated observables 
at hadronic machines include an extra term proportional to the coefficients h % or s*. This 
is similar to analogous relations derived for B s —y <j>{—> K + K~)n + n~ 1301, B s —y VV [32 
and B s —> /i + p _ 31 . The corresponding expressions for time-integrated observables are 


given in Eqs. (42)-(45). However, it seems difficult to extract hi or Sj using time-integrated 
observables, as they are suppressed by small meson-mixing parameters. 

On the other hand, a time-dependent angular analysis with flavour tagging paves the 
way for the observables s*. We identify s 8 and sg as the most interesting observables, as 
they are expected to be large even in the absence of CP violation. We have demonstrated 
that these observables contain new information compared to the angular coefficients Jj, 
and we have built “optimised” versions of these observables with reduced sensitivity to 
form factors. We have focused on two such observables, called Qg and Qg and defined 
in Eqs. (|56]),(57). These observables can be predicted in the SM with good precision (see 
Fig. [Tj) , and show good sensitivity to particular New Physics scenarios (see Figs. [2] and [3]). 

Current analyses of b —> sii transitions point towards deviations compared to SM 
expectations, explained via large NP contributions to Cg (and potentially smaller con¬ 
tributions to other Wilson coefficients). It is particularly interesting and useful to cross 
check this trend from other sources. Our analysis shows that additional information could 
come from time-dependent angular analyses of tagged B ds —> V (—* M 1 M 2 )ff decays, with 
M 1 M 2 a CP eigenstate. We thus encourage exploratory studies to determine the experi¬ 
mental feasibility of such analyses, in particular in the context of a high-luminosity flavour 
factory such as Bclle-II. 
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A CP-conjugate kinematics from invariants 

The kinematics of the four-body decay B —> V[—> M 1 (p 1 )M 2 (p 2 )K + (p+K - (P-) ! is com¬ 
pletely specified by four invariant masses (e.g. s+_ = 2p + -p_, si + = 2 p\-p +1 si_ = 2pi-p_ 
and s 2+ = 2 p - 2 ■ p + ), and the sign of ei 2+ _ = efivpxPiP^P+P-, which defines the parity of 
the final state |44|. We set these kinematic invariants to some fixed values s, Si, s 2 , S3, cr 
such that: 

s+_ = s, si + = si, si_ = s 2 , s 2+ = s 3 , sgn(ei 2+ _) = a . (76) 

The kinematics of the CP-conjugated decay B —> V[—> Mi(pi)M 2 (p 2 )]£ + (p+)£ - (p_) 
is specified analogously. Under the condition of CP-conservation, the differential rates of 
both decays (dT and dT) must be equal at CP-conjugated kinematic points: 

dr(s+_ = s, si+ = si, si_ = s 2 , s 2+ = s 3 , sgn(ei 2+ _) = a) (77) 

= df(s+_ = s, sj_ = si,si + = s 2 , s 2 __ = s 3 , sgn(ei 2 _ + ) = -a) 

= dT (s+_ = s, sj_|_ = s 2 ,Si_ = si,s 2+ = m 2 B - rriy - s - S! - s 2 - s 3 , sgn(e l2+ _) = a) , 

where we have made the replacements {1,2,+,—,} —>• {1,2,—,+} and a —> —a to ac¬ 
count for C and P transformations respectively. In addition, we have used momentum 
conservation and neglected light-meson and lepton masses to write 

s 2 _ = m| - rriy - s+_ - si + - si_ - s 2+ . (78) 


The angular distribution is obtained by expressing these rates in terms of s + _, two 
polar angles 9m, 9f, (0 < 9i < 7r), and one azimuthal anglcj^] 0 (0 < </> < 27 t). In the case of 


dT, the angles are usually defined as 45 


$m = 9\\ Polar angle between the momenta p B and p\ in the rest-frame of V. I 11 
terms of momentum invariants, we find: 


cosdi = 


m B — rriy — s + — 2si + — 2si 


y/(rn 


B 


m 


v 


s_|_) 2 — 4myS + _ 


= ce(s + _,si + ,si_) , 


(79) 


and sin#! = +y/1 — cos6 \ by definition. 


= 9 + : Polar angle between the momenta p B and p + in the dilepton rest-frame. 
In terms of momentum invariants, we find: 


cosd + = 


m 


B 


m 


v 


s + — 2 s 


1 + 


2 s 


2 + 


\/ ( m% - m 2 v - s + _) 2 - 4m|s + _ 
and again sin 9 + = + a/ 1 — cos 6 \. 


= Cg(s + -, Si+, S 2+ ) , 


(80) 


The kinematic angle (f) should not be confused with the mixing angle, not appearing in this appendix. 
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• 0 = 0i2+~: Oriented angle between the planes specified by (pi,p 2 ) and (p + ,/ 7 _ ) in 
the 5 -meson rest frame. The orientation is specified by (in this frame): 

(pi x p 2 ) • (p+ x p_) (p+ + P-) ' [(pi x P2) x (p+ x p_)\ 

COS 0 = ----pz-, Sill 0 = -—-—-—-- — -—-zn- , ( 81 ) 

|Pi X P2I ■ \P+ X p_| \p+ +p-\ • |pi x p 2 | ■ |p+ x p_| 

which in terms of momentum invariants gives: 

cos0i2 + _ = c 0 (s,si + ,si_,s2 + ) , sin012-1— = sgn(ei2+_)y/l - c\ , ( 82 ) 

where c^(s, si+, Si_, s 2 +) = a/( 2 mybc) with 

a = rriy{s + s 1+ ) - (■ m% - s ) [s 1+ (-m 2 B + s + s x _ + s i+ ) + (si_ + si+)s 2+ ] + 
my [s 2 - 2s 1+ + (si_ + si+)(si+ + s 2 _|_) - s (m| - 2si_ - 2 s 1+ - 2s 2+ )] , 


b — y si_ + Si+ — [priy + Si_ + si+)(s + S\- + Si+) , 

c = sjs [(si+ + s 2 + ){rn 2 B - s - s 1+ - s 2+ ) - m 2 v (s + s 1+ + s 2 +)] 
We note that sin0 is proportional to sign(ei 2 +_). 


(83) 


Other possibilities are d M = d 2 or 6 £ = 9_, obtained from (79), (80) by obvious 
replacements. In the case of the CP-conjugate mode, the angles 5 , 9±, 0^ ±T are defined 
analogously. 


In terms of the angular distribution, the CP correspondence in Eq. (77) depends on 


how the angles are defined for dr, relative to dr. We recall that in the case of flavour- 
specific (“self-tagging”) modes, such as B d —» K*(—> K + tt~)££, one might choose any 
convention for dr and df independently, as the final states are different and distinguish¬ 
able. However, this is not the case for untagged flavour-non-specific decays, where the 
final states arising from the B and the B decay cannot be distinguished. We consider 
three different conventions: 


A. dr(s,0 1 ,0 + ,0 12+ _) and df(s, 9\, 9 + , 4>i 2+ _): This is the usual theory convention in 
B d — > K*££, where in both CP-conjugated modes 9m is defined with respect to 
the kaon, 9i with respect to the positively-charged lepton, and the orientation of 
0 is given by 0^ 7r+ _. This is also the only possible convention in untagged decays 
with Mi = Mi and M 2 = M 2 , such as B d — > K*(—t Ks^°)£ + £~ and B s — » 0(—> 


KsKl)£ + £ at hadronic machines. With this convention, Eq. (77) implies 


dr(s+_ = s, cos 61 = Ce(s , si, s 2 ), cos 6 + = c e (s , si, s 3 ), 

cos0 = C0(s,si,s 2 ,s 3 ),sgn(sin0) = a) = 


(84) 
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dT(s + _ = s, cos 6*1 = c e (s, s 2 , s i), cos0 + = c e (s , s 2 ,m| - rriy - s - s\ — s 2 — s 3 ), 
cos 0 = c^(s, s 2 , Si,m 2 B - rriy - s — s x - s 2 — s 3 ), sgn(sin 0) = a) 

We note the following relations: 

Ce{s, Si, s 2 ) = X = cg(s, s 2 , Si) , (85) 

C0(s,si,s 3 ) = Y = -cg(s,s 2 ,m 2 B -rriy- s- si~s 2 - s 3 ) , (86) 

c<fr(s, Si, s 2 , s 3 ) = Z = -c^(s, s 2 , Si, m% - m 2 v - s - s L - s 2 - s 3 ) . (87) 

Therefore, 

rfr(s + _ = s, cos = X, cos 0 + = Y, cos (j) = Z, sgn(sin 0) = er) = 

dT(s + _ = s, cos 0i = X, cos 0 + = — Y, cos 0 = — Z, sgn(sin 0) = a) . (88) 

With the angles defined in this way, the two angular distributions are written as: 

dr = Y^ 9m i 0) * dT = J2 Om, 0) , (89) 

i i 

with Cl,2,3,4,7 — 1 > C-5,6,8,9 = — 1- 

B. dT(s, 0i, 0+, 012+-) and dT (s, 0p 0_, 0i 2+ _): This is the usual experimental conven¬ 
tion for B d —> K*M, where in both modes 9m = 9k and 0 = 0^ 7r +_, but for 9i one 
takes £ + or t~ for the B and B decay respectively. We have: 

dr (s + _ = S, COS 01 = Cg(s, Si,S 2 ), COS 0 + = Cg(s, Si, S3), (90) 

COS0 = c^(s,si,s 2 ,s 3 ),sgn(sin0) = a) = 

dr (s + _ = S, COS 01 = Cg(s, S 2 , Si), COS 0_ = Cg(s, Si, S3), 

COS 0 = c^(s,Si,s 2 , S3), Sgn(sin0) = a) 
which means that, with this convention, 

dr = Ji(s)fi(9 e , 9 m , 0) , dr = Y, Ji(s)fi(9 e , 9 m , 0) • (91) 

i i 

C. dT{s, 0i, 0+, 0i24—) and dT (s, 0 2 , 0 + , 0 2 j + _): This is the only possible convention for 
the case of untagged decays where Mi = M 2 and M 2 = Mi, as for example the 
decay B s —» 0(—> K + K~)i + i~ at a hadronic machine. In this case, 

0T(s_|__ = S, COS 01 = Cg(s, Si, S 2 ), COS 0+ = C e (s,Si,S 3 ), 

COS0 = c^(s,si,s 2 ,s 3 ),sgn(sin0) = a) = 
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( 92 ) 


dT(s+- = s, cos 0 2 = cg(s, m 2 B — rriy — s — si — s 2 — s 3 , S 3 ), 

cos 9 + = c e (s, s 2 , m| - rriy - s — s 1 — s 2 — s 3 ), 
cos 0 = c^(s, m 2 B — my — s — si-s 2 ~ s 3 , s 3 , s 2 ), sgn(sin 0) = -a) . 
Using the relations 

Cfl(s, Si, s 2 ) = X = -c<?(s, m| - my - s - Si - s 2 - s 3 , s 3 ) , (93) 

c e (s, Si, S 3 ) = Y = -Cg(s, s 2 , m| - my - s - Si - s 2 - s 3 ) , (94) 

c</,(s, si, s 2 , s 3 ) = Z = +C 0 (s, m| - m(, - s - si - s 2 - s 3 , s 3 , s 2 ) , (95) 

we have 

dr(s + _ = s, cos 6 *i = X, cos 6 + = Y, cos <j> — Z, sgn(sin 0) = a) = 


dr(s + _ = s, cos $2 — — X, cos 9y = — Y, cos0 = Z, sgn(sin0) = — a) . (96) 
With this convention, the differential rates are given by: 

dT = Ms)fi(9t, 9m, 0) , df=J2 CiMs)fi(0e, 9m, 0) , (97) 


with Ci, 2 , 3 , 4,7 — 1) C5,6,8,9 = — 1- This is the same as Eq. (89) but for a different 
reason. 


We see that conventions A and C yield the same relation between dr and dT, but they 
apply to different kinds of modes. For decays into flavour-specific modes, convention B is 
also possible, but with a different relationship between dr and dT. In the present paper, 
we choose convention A for decays into flavour-specific modes as well as for B —> V(—> 
M\MT)ti decays with M\ = Mi, M 2 = M 2 , and convention C for B —> V(—> MiM 2 )££ 
decays with Mi = M 2 , M 2 = M\. 


B CP-parities associated to transversity amplitudes 

We consider the decay B —>• VN (cf. Eq (J4]) ) where V and N are unstable particles, V 
decaying into two particles M\ and M 2 . As shown in Ref. |35], the CP-parity of a final 
state X is given by 

Vx = C(-l r (98) 

with r = r(Mi) + r(M 2 ) + t(N) the “transversity” of the state M 1 M 2 N (defined below), 
and C depends on the class of decay: 

• class 1: V (not necessarily with a definite spin) decays into Mi and M 2 which are 
CP-eigenstates, and N decays into a CP-eigenstate. In this case, 

£ = t)(N)t](Mi)t](M 2 ). (99) 
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X 

N 

V (N) 

S 

t(N) 

Vx 

0 

vector 7 *, Z (or axial) with eo 

1 

1 

1 

V 

II 

vector 7 *, Z (or axial) with e|| 

1 

1 

1 

V 

_L 

vector 7 , Z (or axial) with ej_ 

1 

1 

0 

-77 

5 

scalar H 

1 

1 

0 

-v 

t 

vector 7 *, Z (or axial) with e t 

-1 

0 

0 

V 

t 

pseudoscalar A 

-1 

0 

0 

77 


Table 3: Properties of the transversity amplitudes involved in B —>■ V££: CP-parity of N, 
spin of the lepton pair, transversity of N, and CP-parity of the final state. 

• class 2 : V (with a definite spin sy) decays into spin-0 Mi and M 2 which are CP- 
conjugates, and N decays into a CP-eigenstate. In this case, 

i = v {N){-iyy. ( 100 ) 

• class 3: V and N are CP-conjugates with a definite spin sy, with V decaying into 
spin-0 Mi and M 2 . In this case, 


£=(-!)*''• 


( 101 ) 


Here 77 (H) and sy are the intrinsic CP-parity and spin of the particle V. The first class 
is illustrated by the time-dependent analysis of Bd —> Ks'k °) [361. For the 

class-1 processes ( Bd —» K*(—> Ksx°)££ and B s —> <j >{—* KsKl)££) and class-2 process 
(B s —» </>(—>■ K + K~)££) of interest, we have 

Vx = ?/(iV)(-l ) T(Ar)+1 77 (102) 


where r) = — rj(Mi)rj(M 2 ) for class-1, and 77 = 1 for class-2. For all the processes considered 
here, the combinations of intrinsic CP-parities yield 77 = 1 . 

In order to determine the CP-parity of the different transversity states, we have thus 
to determine r](N) and r(N): 


Using the language of Ref. 35 , we see that Aq, A±, A\\,As t t are respectively associ¬ 
ated with the combinations of hclicity amplitudes denoted 0 , Q\ 0 o> @10 0 @0 oo> 
with respective CP parities —£, —For our decays, it implies that we should 
have the following associations (modulo 2 ): 


A), A\\ \ t = 1 , A ± ,A tl A s : T = 0 (103) 

The states corresponding to different transversities can be accessed through the 
angular analysis of the decay. 
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• r](N ) can be determined from the assumed quantum numbers of the intermediate 
boson (spin, polarisation, parity). rj(N) is identical for spin-1 particles with vector 
or axial couplings, whereas scalar and pseudoscalar N have opposite CP-parities. 
One can also notice that this constrains the spin of the emitted £ + £~ pair. If we 
denote its angular momentum l and total spin s (either 0 or 1), the P-parity of such 
a fermion-antifermion pair is given by (—1) /+1 , its C parity by (—l) z+s , so that its 
CP-parity is (—l) s+1 . Since we assumed that the decay N —» £ + £~ conserves CP- 
parity, we have rj(N) = (—l) s+1 . We have l — s — 0 corresponding to a pseudoscalar 
N, l = s = 1 corresponding to a scalar N, l — 0, s — 1 corresponding to a (real) 
vector/axial N, l — l,s = 0 corresponding to a time-like vector/axial N (this can 
be checked from the CP-parity of the corresponding fermion-antifermion currents). 

For each amplitude, we can determine the intermediate virtual boson N with the 
appropriate quantum numbers, the corresponding spin of the lepton pair, the transversity 
associated, and the CP-parity of the final state, as indicated in Table [3j We see in 
particular that we agree with the assignments for the class-1 decay Bd —» J/iJjK*(— y 
/ls7T°) 136 . In the end, we have 

r]x = V for X — LQ, L\\, RQ, R\\,t ; rjx = —r/ for X = L _L, R _L, S . (104) 

We impose t(N) = 0 for a vector N with timclike polarisation, to obtain the same CP- 
parity as in the pseudoscalar case. This agrees with the expectation that A t should have 
the same CP-parity as Aq. 


C Expressions for the coefficients Si and hi 

The coefficients s t are given by 


s u = 2 ^Y^lm[e i<l, {A I [A I [* + Aft Aft* + AfAf + Aft Aft*}} (105) 

1 rn ^ 

+ -^Im ^{A^Af* + Aft Aft*} - {Aft Aft* + AftAft*}\ 

sic = 21m [e i(t> {AftAft* + Aft Aft*}) (106) 


S2s ~ 

S2c — 
S 3 = 


8m 

q2 


2 r 


Im [e i(t> {A t A* t }\ + Im[e^^ - Aft Aft*) + 2/3 f 2 Im[e^kl s klI 


R2 

^Im [e**{AftAft* + Aft Aft* + Aft Aft* + Aft Aft* }] 

-2/^Im [e^{AftAft* + AftAft*}) 
ffilm[e irp {AftAft* - Aft Aft* + Aft Aft* - 
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(107) 

(108) 
(109) 







54 = + A * A \t *> - ^{AftAft* + Aft Aft *}] 


«5 


= 


Im[e # {^Af - Aft Aft*} - e~ i4, {AftAft* - Aft Aft*}] 


mg 


Im [e«{AftA* s + AftA* 8 } - e^*{AftA* s + 


s 6 , = 2^Im[e^{Af Ai* - Aft Aft*} - {Aft Aft* - Af Af}] 
S6c = 4/^Im[e^{l^ + - e^K^s + 4 }%,}] 


s 7 — —V2/3e 


Re[e i4> {AftAft* - Aft Aft*} - - AM*}] 




Re[e^{^^ + AftA* s } - e~* {AftA* s + 4^}] 


-=/9flte[e^{A^i* + }] 

v2 


*8 = 

s 9 = /3jFte[e^{^84i' + Af } - e^A^A? + Affylf}] 


( 110 ) 

( 111 ) 


( 112 ) 

(113) 

(114) 


(115) 

(116) 


The coefficients h t are given by 

h ls = —^Re[e l4> {A L L A L ^ + Aft Aft* + AftAft* + Aft Aft* }] (117) 

1777 ^ ^ 

+— ^Re[e i(t> {AftAft* + AftAft*} + {AftAft* + Aft Aft*}] 

h\ c = 2Re[e i< t > {Aft Aft* + A* A**}] (118) 

+^[Re[e^l^*] + ^{e** Aft Aft* + Aft Aft *}] + 2$Re[e^{l s ^}] 

«2 __ _ __ __ 

/i 2s = ^Re[e^{7L^* + Ajj'Aji'* + AftAft* + Aft Aft* }] (119) 

h 2c = -2f3 2 R ^{AftAft* + A* Aft*}] (120) 

h 3 = /3jR,e[e^{AftAft* - Aft Aft* + AftAft* - Aft Aft* }] (121) 

h A = ~^= p 2 Re[e ir ^ {Aft Aft* + Aft Aft*} + {Aft Aft* + Aft Aft* }] (122) 

v 2 
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— y/2,fie 


Re[e^{Af Af* - A*Af} + e 


^{Aq Af* 


— a r a r * 

/1 0 /lj_ 


}] 


h($s — 
hec — 


-^LRe\e‘*{A^A‘ s + XfiUj} + e^A^ + ^J}] 

2&Re[e^{Af Af* - If Af } + e^jAf Af* - Af Af }] 
4A^ 5 Re[e‘»{i 0 M* s + l«,4* s } + e-*K4 + <lj}] 


/17 — V2(3e 


Im[e* 




a r a r * 
^0 Ai 


} + 


*^{Af A L * 


AfAf* 


}] 


+^=lm[e i4, {A^_A* s + A±A* S } + e~ ill> {A L L A* s + Af A^}] 

/i 8 = 4=$Im [e^{Af Af* + If Af } + e^Af Af + Af Af }] 

V 2 

ft, = + Ifvlf } + + j4f j4f*}] 


(123) 


(124) 

(125) 

(126) 


(127) 

(128) 


In the above expressions, the amplitude Ax denotes the amplitude Ax{B —1 /), 
without applying CP-conjugation to the final state. One has the relation 

A x = VxAx (129) 

where A x can be obtained from A x by changing the sign of all weak phases. 
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